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Abstract.  We  continue  the  study  on  equilibria  of  the  Smoluchowski  equation 
for  dilute  solutions  of  rigid  extended  (dipolar)  nematics  and  dispersions  under 
an  imposed  electric  or  magnetic  field  na.  We  first  provide  an  alternative 
proof  for  the  theorem  that  all  equilibria  are  dipolar  with  the  polarity  vector 
parallel  to  the  external  field  direction  if  the  strength  of  the  permanent  dipole 
(fi)  is  larger  than  or  equal  to  the  product  of  the  external  field  (E)  and  the 
anisotropy  parameter  (qq)  ^  ^  |q^o|F).  Then,  we  show  that  when  < 

|q:o|F,  there  is  a  critical  value  o*  >  1  for  the  intermolecular  dipole-dipole 
interaction  strength  (a)  such  that  all  equilibria  are  either  isotropic  or  parallel 
to  the  external  field  if  q  <  o*;  but  nonparallel  dipolar  equilibria  emerge  when 
a  >  a* .  The  nonparallel  equilibria  are  analyzed  and  the  asymptotic  behavior  of 
a*  is  studied.  Finally,  the  asymptotic  results  are  validated  by  direct  numerical 
simulations. 


1.  Introduction.  Extended  rigid  nematic  polymers  are  consisted  of  dipolar 
rigid  nematic  polymers  which  carry  electric  dipole  or  magnetic  moments  perma¬ 
nently  or  under  an  imposed  electric  or  magnetic  field  [siiniEniiii  Particle  dis¬ 
persions  of  electric  or  magnetic  moments  in  solutions  are  important  materials  that 
have  many  applications  in  industries  as  well  m  For  example,  in  the  manufacture 
of  flexible  magnetic  data  storage  media,  a  fast-moving  substrate  is  coated  with  one 
or  more  layers  of  dispersions  (or  inks)  composed  of  sub-micrometer-sized  magnetic 
particles.  These  magnetic  particles  are  coated  with  a  polymeric  binder  and  dis¬ 
persed  in  a  solvent.  To  achieve  high  data  storage  capacity,  it  is  very  important  to 
control  particle  orientation  in  the  inks  which  is  sensitive  to  external  aligning  fields 
(flow  and  magnetic). 
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Doi  and  Edwards  used  the  Doi-Hess  kinetic  theory  to  model  the  homogeneous 
flows  of  semi-dilute  rodlike  nematic  polymers  that  may  possess  electric  dipoles  [3 
cni  nni-  Recently  the  Doi-Hess  kinetic  theory  was  extended  to  model  magnetic 
dispersions  in  viscous  solvent  PEI  EH  El.  Both  theories  are  essentially  identical. 
Without  loss  of  generality,  we  focus  on  the  dipolar  rigid  nematics  and  magnetic 
dispersions  that  are  sensitive  to  magnetic  fields.  Using  a  closure  approximation 
a  mesoscale  constitutive  model  for  magnetic  dispersions  was  developed  to  study 
material  functions  for  shear  flows  and  phase  transitions  at  equilibrium  P  P .  The 
phase  transitions  were  re-visited  by  solving  the  Smoluchowski  equation  exactly  in 
equilibrium  In  this  paper  we  extend  our  earlier  results  in  OP  ESI  CD  to  study 
equilibrium  solutions  of  the  Smoluchowski  equation,  in  which  the  the  first  moment 
vector  (polarity  vector)  is  not  parallel  to  the  external  field.  We  will  provide  a 
detailed  study  on  properties  of  the  non-parallel  solutions. 

The  paper  is  organized  as  follows.  In  Section  2  we  give  a  brief  description  of  the 
properties  of  magnetic  dispersions  under  an  external  magnetic  field.  After  reviewing 
the  equilibrium  solution  of  Smoluchowski  equation,  we  prove  a  series  of  theorems 
that  characterize  the  properties  of  the  equilibrium  solutions  whose  polarity  vector  is 
not  parallel  to  the  external  field.  The  asymptotic  result  for  the  critical  value  of  the 
dipole-dipole  interaction  strength  is  derived  in  Section  3.  The  details  are  put  in  the 
Appendix.  Numerical  results  are  presented  in  Section  4  to  validate  the  asymptotic 
results.  Finally,  in  Section  5  we  summarize  the  main  results  of  the  paper. 

2.  Equilibria  of  Smoluchowski  equation  for  magnetic  dispersions  or  rigid 
nematic  dipolar  polymers  under  an  imposed  magnetic  field.  The  ori¬ 
entational  probability  density  function  for  homogeneous  flows  of  magnetic  nematic 
dispersions  or  rigid  nematic  polymers  p(m,  t)  is  governed  by  the  Smoluchowski 
equation  ESI  It  has  been  shown  that  the  equilibrium  solution  of  the  Smoluchowski 
equation  is  given  by  the  Boltzmann  distribution  0 13  0  El  EH  Ea  ESI  mi  ESI 

p(m)  =  ^exp[— [/(m)],  Z  =  /  exp[— t/(m)]dm  (I) 

Z  Js 

where  m  represents  the  permanent  dipole  direction  of  the  nematic  molecule  in 
dispersions  when  /i  ^  0  or  the  induced  dipole  direction  when  /i  =  0,  U(m)  is  the 
total  potential  energy  normalized  with  respect  to  kT  {k  is  the  Boltzmann  constant 
and  T  the  absolute  temperature),  Z  is  the  normalizing  constant  or  the  total  partition 
function,  and  S  denotes  the  unit  sphere  S  =  {m|||m||  =  I}. 

For  dilute  solutions  of  dipolar  rigid  nematic  polymers  under  an  imposed  magnetic 
field  (E),  the  total  potential  energy  density  is  given  by 

Qin 

t/(m)  =  — a(m)  •  m  —  fiE  ■  m - (2) 

where  uq  is  the  difference  of  the  polarizability  parallel  and  perpendicular  to  m 
known  as  the  anisotropy,  fi  is  the  strength  of  the  permanent  dipole  moment,  a  is 
the  strength  of  the  intermolecular  dipole-dipole  interaction,  (m)  is  the  first  moment 
of  the  pdf  also  known  as  the  polarity  vector, 

(m)  =  /  nip(m)dni.  (3) 

Js 

In  this  paper,  we  study  dilute  nematic  polymer  solutions  or  dispersions  so  the 
excluded  volume  effect  is  neglected  in  the  potential  energy.  Given  the  material 
parameters  (a,  ag,  /i,  and  E),  the  potential  energy  as  well  as  the  equilibrium  solution 
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are  completely  determined  by  (m)  provided  the  non-linear  integral  equation  for  (m) 
©  is  solved.  We  shall  focus  our  attention  on  the  solutions  of  equation  lEJ  for  various 
parameter  regimes  in  this  study. 

We  first  establish  a  coordinate  system  for  m;  we  select  the  z-axis  to  be  the 
direction  of  E,  the  y-axis  to  be  perpendicular  to  the  plane  spanned  by  (m)  and  E; 
we  also  select  the  positive  directions  of  the  x-axis  and  y-axis  such  that  (m)  has  a 
non-negative  x-coordinate  (if  (m)  and  E  are  parallel  to  each  other,  then  the  choice 
for  the  x-axis  and  y-axis  is  not  unique).  It  is  noteworthy  that  for  a  given  coordiate 
system  with  the  z-axis  parallel  to  the  external  field,  due  to  axisymmetry  of  the 
problem,  if  (ri,0,  ra)  is  a  solution,  then  (ri  cos6>,  ri  sin0,  xa)  is  also  a  solution  for 
any  value  of  0.  Here  we  select  (ri,0,  ra)  with  xi  >  0  to  represent  this  group  of 
solutions. 

In  this  Cartesian  coordinate  system,  we  have 

m  =  (mi,TO2,ma),  E  =  E(0,  0,1),  (m)  =  (xi,  0,  xa),  xi  >  0, 

ao 

t/(m)  =  — axi  TOi  —  {fiE  +  axa)ma - ^  ^3" 

The  nonlinear  integral  equation  ©  becomes 

xi  =  /  mi/9(m)dm, 

Js 


0  = 


m2p(m)(im. 


(5) 


Xa  =  /  ma/9(m)dm. 

Js 

Note  that  our  choice  of  the  Cartesian  coordinate  system  makes  E  positive.  In  this 
paper  we  consider  the  case  where  a  >  0,  E  >  0,  |ao|  >0  and  fj,  >  0.  In  [2S1  we 
have  shown  that  if  a  <  1,  then  (m)  must  be  parallel  to  E.  Furthermore,  we  have 
shown  that  if  /r  >  |ao|£';  then  (m)  must  be  parallel  to  E.  In  this  paper,  we  will 
first  give  an  alternative  proof  for  this  result  based  on  the  torque  balance  argument 
and  then  continue  on  analyzing  the  behavior  of  solutions  whose  polarity  vector  is 
not  parallel  to  the  external  field  direction. 

For  the  theorems  below,  we  introduce  a  spherical  coordinate  system  in  which  the 
y-axis  points  to  the  north  pole.  We  use  {ip,  ()  to  denote  this  spherical  coordinate 
system  where  ip  is  the  polar  angle  and  ({  is  the  azimuthal  angle.  Later  on  in  this 
paper,  we  will  use  another  coordinate,  which  we  denote  by  {(p,  9)  for  distinction,  in 
which  the  z-axis  points  to  the  north  pole. 

In  the  {ip,C)  spherical  coordinate  system,  we  have 


{mi,m2,m^)  =  (sini/isinC,cosi/’,sin'0cosC), 

u{iP,0  =  UmuM,0  +  Ue.M,0, 


(6) 


where  the  mutual  interaction  and  the  external  part  of  the  potential  are  given  re¬ 
spectively  by 

UMuti'P,  C)  =  xi  sin')/; sin C,  —  ar^  sin')/; cos (p, 


UExt{'ip,C) 


— /r  E  sin  ip  cos  (p  —  ^  E^  sin^  ip  cos^  p. 


(7) 


910 


HONG  ZHOU,  HONGYUN  WANG  AND  QI  WANG 


The  equilibrium  pdf  is  given  by 
1 


Z  = 

Jo  J--K 

The  nonlinear  integral  equation  for  ri  and  becomes 


exp[— 17(V',C)]^Csin'!/'di/’-  (8) 


=  Jo 


'^3  =  JJ  11^  smip  cos  Cp{ip,Odil}smil}dC. 

Theorem  1.  At  equilihrium,  the  external  torque  about  the  y-axis  vanishes,  i.e., 

d 


(9) 


9C 


=  0- 


Proof.  We  first  show  that 


)  =  0. 


(10) 


(11) 


Physically,  this  quantity  is  the  (negative)  total  torque  on  the  system  about  the 
y-axis.  Since  the  system  is  in  equilibrium,  the  total  torque  should  be  zero,  i.e., 


lo  J  - 
Next  we  have  that 
d 


1  r  r  d 

'Z  j  J  -  tliSxt('*/',C)] '^Csini/'dV' =  0. 


(12) 


C)  )  =  (a^3  sini/jsinC  -  ori  sini/jcosC)  =  ar^ri  -  arir^  =  0,  (13) 

dC  / 

where  we  have  used  ®.  Equations  113  and  m  lead  immediately  to  m-  We  note 
that  m  is  a  consequence  of  the  Newton’s  third  law  for  mutual  interaction.  □ 


Theorem  2.  If  p>  laolE,  then  the  polarity  vector  (m)  of  an  equilibrium  solution 
is  non-zero  and  must  be  parallel  to  E. 


Proof.  We  need  to  prove  that  all  solutions  of  satisfy  ri  =  0  first.  We  prove  it 
by  contradiction.  Suppose  there  is  a  solution  of  satisfying  ri  >  0  (recall  that 
the  coordinate  system  is  selected  to  make  ri  non- negative) .  We  are  going  to  show 

that  p  >  \ao\E  and  ri  >  0  lead  to 

cm  proved  in  Theorem  1  above. 

Exploiting  the  fact  that  UExti'4’,  C)  is  an  even  function  of  f,  we  have 


A 


UExtiifX)  )  ^  Oj  which  contradicts  with 


(^t/u.t(V’,C)) 

^  ^  Jo  J  ^^^“*(^’‘^Hexp[-f7(t!’,C)]  -  exp[-t7(V’,-C)]}rfCsin'0d'0. 

In  the  above,  the  first  factor  of  the  integrand  satisfies 


d 

-^UExtiipX)  =  pEsintpsinf  -|- sin^  V' cos  (sinC 


(14) 


(15) 


=  (/r -I- aoEsini/jcosC-E'Sin'f/'sin^  >  0  forCG(0,7r),  '0G(O,7r). 
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The  second  factor  of  the  integrand  satisfies 
exp[-C/(V',  C)]  -  exp[-C/(V',  -()] 


sm‘‘ Ip  cos^  C  +  {fiE  +  ar 3)  sin'll)  cos  sinh(ari  sini/jsinC)  (16) 


=  2exp 

>0  for  C  G  (0,  tt),  Ip  G  (0,  tt). 
Substituting  03)  and  into  m  we  obtain 


A 


UExti'P,  C)  /  >  Oj  which  contra- 


□ 


diets  with  m-  Therefore,  when  >  |ao|i?,  (m)  must  be  parallel  to  E. 

Next,  we  have  to  show  that  the  polarity  vector  is  not  a  zero  vector,  i.e.,  7^  0. 

It  follows  from  0  that 

ra  =  (sin '0  cos 

(17) 

=  I  sin^  ip  cos  C  sinh((ar3  +  ^E)  sin  1/;  cos  C)e^^'  sin^  v-cos^  ‘^dp^dC. 

This  implies  that  ^  0  since  ra  =  0  contradicts  to  the  above  equation. 


Now  let’s  look  at  the  intuitive  physical  picture  of  the  non-parallel  solution.  Each 
polymer  is  subject  to  two  potentials:  i)  the  external  potential  caused  by  the  mag¬ 
netic  field  and  ii)  the  mutual  interaction  potential  caused  by  the  mean  field.  The 
net  torque  effect  of  the  mutual  interaction  potential  is  zero  by  Newton’s  third  law. 
Therefore,  it’s  the  torque  due  to  external  field  that  affects  the  mesoscale  orientation. 
In  addition,  it  has  been  shown  that  a  dipolar  nematic  phase  forms  beyond  =  3 
when  the  material  is  absent  of  any  external  field.  In  the  spherical  coordinate  system 
where  the  pole  is  on  the  z-axis,  the  external  potential  is  given  by 

UExt{(p,d)  =  - fiE  cos  (p-  ^E^cos^(/)  =  -oqE'^  ’ 


which  has  two  stationary  points  (p  =  0  and  </>  =  tt  when  fi  >  |ao|E.  When  fj,  <  |q!o|E, 
however,  the  external  potential  has  a  third  stationary  point,  (po,  given  by 


cos{(po) 


d 

aoE 


The  torque  due  to  external  potential  vanishes  at  the  stationary  point,  leaving  the 
possibility  for  a  nonparallel  equilibrium  at  /r  <  |ao|i?  and  a  large  enough. 

A  necessary  condition  for  the  existence  of  the  non-parallel  solution  indicated  in 
Theorem  2  above,  <  |ao|E,  seems  to  suggest  that  for  large  external  field  (E),  the 
non-parallel  solution  shall  always  exist.  However,  the  theorem  below  (Theorem  3) 
tells  us  that  when  oq  >  0,  non-parallel  solution  does  not  exist  for  large  external 
field.  This  is  physically  reasonable.  When  oq  >  0,  the  effect  of  the  external  field 
both  directly  on  permanent  dipole  and  indirectly  on  the  polarized  dipole  tends  to 
align  nematic  polymers  with  the  external  field  (parallel  and  anti-parallel).  As  one 
will  see.  Theorem  3  below  does  not  exclude  the  existence  of  non-parallel  solution 
for  large  external  field  when  ag  <  0.  Our  numerical  simulations  show  that  when 
ao  <  0,  non-parallel  solution  alway  exists  for  large  external  held.  When  ag  <  0, 
the  effect  of  the  external  held  on  the  polarized  dipole  tends  to  align  polymer  rods 
perpendicular  to  the  external  held,  which  is  what  we  observe  in  the  numerical  sim¬ 
ulations. 
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Theorem  3.  If  (ri,r3)  is  a  solution  of  with  ri  >  0,  then  +  ar^)  <  0. 

Proof.  We  prove  the  theorem  by  contradiction.  Suppose  there  is  a  solution  of  @ 
satisfying  ri  >  0  and  aoifJ-E  +  ar^)  >  0.  We  will  show  that  ri  >  0  and  aoifrE  + 

ars)  >  0  lead  to  (^■^UExt{'f’,o'^  >  0,  which  contradicts  with  (tinil  in  Theorem  1. 
First  we  have 
d 

tttUexM,  C)  =  /ri?  sin i/j  sin  ^  +  aoE^  sify  if  cos  C  sin 

dC 

>  aoE"^  sin^  if  cos  C  sin  C  for  C  €  (0,  tt)  ,  tp  £  (0,  tt)  . 
Substituting  it  into  d  and  using  d,  we  obtain 

{§-^v,M:0)  >  fy 

pTT  nTZ 

/  /  aoE^  sin^  cos  C  sin  ^  {exp[— 17(^,  ()]  —  exp[— [/(t/j,  — C)]}  (18) 

^0  ^0 


/o  ^0 

I 


t/2 


4i?^  sin^  ip  cos  C  sin  (  exp 


0  ^0 


^i?^  sin^  ?/;  cos^  C  9  {i’tOdC  sin  ip  dip, 


where 


g{ip,  C)  =  tto  sinh[(/ri?  +  ar^)  sini/>cosC]  sinh(ari  sini/^sinC) 


>  0  for  C  S  (0, 7r/2),  ip  G  (0,  tt),  ao{fJ.E  +  ar^)  >  0. 


Thus,  we  have 


>  0,  which  contradicts  with  d. 


□ 


In  |25|  we  conjectured  and  numerically  confirmed  the  existence  of  non-parallel 
solutions.  In  the  theorem  below  we  give  a  rigorous  proof  for  the  existence  of  non¬ 
parallel  solutions. 

Theorem  4  (Existence  of  non-parallel  solutions).  For  /r  <  |ao|i?,  there  exists  an 
a  (depending  on  g,,  ag,  and  E)  such  that  when  a  >  d  there  is  a  solution  of  I^) 
satisfying  ri  >  0. 

Proof.  Consider  the  free  energy  density 

F(a,r3,ri)  =  log  J  exp[-U{ip,C)]dCsmipdip'^  -  ^{rf  +  rj).  (19) 

Equation  is  equivalent  to  the  Euler-Lagrange  equation  derived  from  the  free 
energy  density: 

— F(a,  r3,  n)  =  0,  r3,ri)  =  0.  (20) 

We  express  and  ri  in  polar  coordinates 


r3  =  r  coscj, 


n  =  r  sincj. 


(21) 


NONPARALLEL  SOLUTIONS  OF  EXTENDED  NEMATIC  POLYMERS 


913 


The  function  F{a,r3,ri)  becomes 
F{a,r,ui) 


log(/  J  exp[arsintpcos{C  -  u;)  -  UExt{ip,0]dCsiniljdijj^  -  ^22) 

log  J  exp[arsin^cos(C)  -  UExt{tp,C  +  uj)]d(  sin 


Equation  is  equivalent  to 
d 


dr 


F{a,r,u;)  =  0, 


duj 


F(a,r,  w)  =  0. 


(23) 


So  we  only  need  to  show  that  for  /r  <  |ao|E,  there  exists  an  a  (depending  on  /r, 
ao  and  E)  such  that  when  a  >  a  there  is  a  solution  of  it^  satisfying  ri  >  0.  We 
prove  it  in  several  steps. 

Step  1  :  Consider  the  function 

^  J^^exp[Xsiml>cosC-UExti'tp,C  +  ^)]dCsimpdip 

Ti{X,io)  = - 


We  have 


Jo  /_^  exp[Asini/i cos C]dC sini/idi/i 


lim  ri(A,w)  =  exp[-UExt{-z,uj)]  =  T2{lu). 

X — *oo  Z 


and  the  convergence  is  uniform  in  lo.  The  uniform  convergence  is  needed  in  Step  5 
below.  We  prove  this  result  using  the  (e,  A)  notations  in  Calculus.  Note  the  e  used 
in  this  step  should  not  be  confused  with  the  e  used  in  the  later  part  of  the  paper. 
For  any  e  >  0,  we  want  to  show  that  there  exists  A  such  that  A  >  A  implies 


|Ti(A,a;)-T2(w)|  <e, 


(24) 


for  all  ta  G  [— tt,  tt].  We  prove  in  several  substeps. 

d  d 

Substep  lA:  We  notice  that  UExt{ip,0,  -^UExtii^X)  and  -^UExtXX)  are 
bounded  for  all  ij}  and  C,.  That  is,  there  exits  M  such  that 


\UExtXX)\<M, 


A 

d-tp 


UexMX) 


<  M, 


A 


UexMX) 


<  M, 


(25) 


for  all  Ip  and  C-  It  follows  that  there  exists  0  <  5  <  ^  such  that  !■!/)—  ^|  <6  and 
Id  <  (5  implies 


TT  1 

lexpi-UExtiipX +  X)  -  <i^I>[-UExt{^,uj)]\  < 


(26) 


for  all  Lu. 

Substep  IB:  Let  us  introduce  two  shorthand  notations  for  the  integration  domains. 

D  = 


{(V’,c)  x-^\<  |,ici  < ttI , 

Ds  =  [i'ipX)  IV’- |l  <  <^JCI  < 
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Result  (EEl  leads  to 

Ids  Asini/’  cos  C  —  UExtii’,  C  +  s\rL'ij)dij) 


fjy  exp[A  sini/j  cos  C]dC  sin'0d'0 


-  exp[-l7Ext(^,uj)] 


1 

<  -e 
4 


(27) 

for  all  Lu. 

Substep  1C:  In  D\Ds,  we  have  either  cosC  <  cos 5  or  sinip  <  cos^  or  both.  It 
follows  that 


lD\Ds 


exp  [A  sin  ip  cos  (]d(  sintpdtp  <  2tt'^  exp  [A  cos  5] 


Thus,  we  have 


lim  Aexp(— A)  /  exp[Asini/>cosC]dCsin'0d'0  =  0 

JD\Ds 


(28) 


(29) 


In  Ds ,  using  the  fact  that  cos  2  >  1  —  letting  =  '/Xi'tp  ~  f )  and  =  a/ACi 


we  obtain 


Aexp(— A)  /  exp[Asin'!/;cosC]dCsinV’di/’ 

J 

'  f )  dC  cos  6dtp 


>  A  /  exp 
Jds 


=  cos  (5 

— >  27r  cos  5  as  A  ^  oo 
Combining  (PI  and  (PI  yields 


P^/xs 

/  exp 

J-^sJ 

-CA5 

2 

dC'diA' 


(30) 

(31) 


Jd\Ds  exp  [A  sin  ■)/)  cos  C]dC  shiipdip 
A^oo  exp  [A  sin  ^  cos  dC,  sin  ipdip 

Substep  ID:  Since  UExt{'4’tO  is  bounded,  we  have 

y  Id\Ds  exp[Asin^cosC  -  UExtiipX  +  uj)]dC sin tp dip 

A^oo  exp  [A  sin  cos  C  —  UExtidP,  C  +  ‘^)]'^C  sinipdip 
fnx  n  exp [ A  sin  ip  cos  C]dC  sinipdip 

<  exp(2M)  hm  j  =  0  (32) 

A^oo  exp[Asin'!y  cosCJuCsini/yai/y 

and  the  convergence  is  uniform  in  lo.  Combining  m  and  1211,  we  obtain  that  there 
exists  A  such  that  A  >  A  implies 

/p  ^  C-c/Ext (b.C+<^)] d(  sin  ipdip 


Jjj  g[Asmi/>cosC](j^  sinipdip 
g[AsinV.cosC-C/E,.t(V’.C+<x^)]rf^sim/,d^ 


(33) 


fos  ^  Ipdip 


1 

<  -e 
4 


for  all  Lu.  Finally,  combining  (HZl  and  1221,  we  conclude  that  A  >  A  implies 

/p  g[A  sin  V-  cos  C-C/Ext  (V'.C+'x^)]  d(;  sin  Ipdip 


Jjj  e  b  cos  c]  sin  ipdip 


-  exp[-UExt{^,i^)] 


1 

“  2^ 


(34) 
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for  all  oj,  which  leads  immediately  to  ll^. 
Step  2:  By  definition 


T2(w)  =  exp[fj,E  cos  uj  +  cos^  lu] 


(35) 


Let  Wo  be  the  angle  satisfying  cos(wo)  = 


aoE 


and  0  <  Wo  <  TT.  T2(w)  attains  a 


minimum  (ap  >  0)  or  a  maximum  (ao  <  0)  at  wq.  Below  we  present  the  proof  for 
the  case  of  ao  >  0.  The  case  of  ao  <  0  can  be  handled  similarly. 

We  choose  <5  >  0  such  that 

T2(wo  —  6)  —  T2(wo)  >  0,  T2{ujo  +  5)  —  T2(wo)  >  0, 


Wo  —  (5  >  0,  ujQ  +  6  <  TT. 

Using  the  result  in  Step  1,  we  have 

lim  [Ti(A,wo  -  S)  -  Ti(A,wo)]  >  0,  lim  [ri(A,  wo  +  (5)  -  ri(A,wo)]  >  0. 

A — ^oo  A — *oo 

Step  3:  From  Step  2,  there  exists  Ao  such  that  A  >  Ao  implies 

Ti(A,  Wo  —  6)  —  Ti(A,  Wo)  >  0,  ri(A,wo  +  5)  —  Ti(A,  wo)  >  0. 

As  a  result,  for  A  >  Ao,  Ti(A,w),  as  a  function  of  w,  attains  a  local  minimum  in 

(wo  -  (5,  Wo  +  S). 

F(a,r,  w)  is  related  to  Ti(A,w)  as 

F(a,r,  w)  =  log(Ti(ar,  w))  +  log  J  exp[ar  sin^  cos  C]dC  sin'i/'d'!/;^  — 

Notice  that  only  the  hrst  term  depends  on  w. 

Step  4:  From  Step  3,  we  see  that  when  a  >  2Ao,  for  each  r  in  [0.5, 1],  F(a,r,  w), 
as  a  function  of  w,  attains  a  local  minimum  in  (wo  —  <5,  wo  +  <5).  When  a  >  2Ao,  the 
location  of  the  local  minimum  defines  a  function  w(r)  for  r  in  [0.5, 1].  w(r)  satisfies 


d 

-—E{a,  r,  w(r))  =  0,  0  <  wo  —  5  <  w(r)  <  wo  +  <5  <  tt. 
aw 


Step  5  :  Consider  the  function 


T3(A,w)  = 


Jo  I-Tv  sin  V' cos  C  exp  [A  sin  ^  cos  C  —  UExt{'4’^  C  +  ^)]^C  sinipdip 
/o’"/.r^exp[Asin^cosC-  UExt{'fp,C  +  w)]dC sin V'# 


(36) 


(37) 


We  have  T3(A,w)  <  1  for  Hnite  A  since  T3(A,w)  can  be  viewed  as  the  average  of 
sin'!/)COS(^  with  respect  to  a  probability  density.  Furthermore, 

TT 

lim  T3(A, w)  =  sin  —  cosO  =  1. 

A — >-oo  2 

Since  UExti’ipjC)  is  bounded,  this  convergence  is  uniform  in  w.  As  a  result,  there 
exists  Ai  such  that  for  A  >  Ai  we  have 

i<T3(A,w)<l,  for  all  w. 

Step  6:  Differentiating  F(a,r,  w)  with  respect  to  r,  we  get 
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Thus,  when  a  >  2Ai,  we  have 

^F(a,  1,0;)  =  a[Tz{a,uj)  -  1]  <  0. 

That  means  when  a  >  2Ai  for  each  lo,  -^F(a,r,ijj)  =  0,  as  an  equation  for  r, 
has  a  solution  in  [0.5,1].  When  a  >  2Ai,  the  solution  defines  a  function  r{uj)  for 
uj  S  [0,  27r].  r(o;)  satisfies 

^F(a,r(w),o;)  =  0,  ^  <  r(w)  <  1.  (38) 

Now  we  combine  TO  and  TO  to  finish  the  proof.  Let  a  =  2max(Ao,  Ai).  When 
a  >  a,  both  (TOl  and  TO  are  valid.  Let  (ra,LJa)  be  the  intersection  point  of  the 
curve  aj(r)  in  (EH)  with  curve  r(aj)  in  EHl),  as  illustrated  in  Figure  1.  At  {va^uja), 
We  have 

d  d 

—F{a,ra,uja)  =  0,  ^F(a,  r^,  o;^)  =  0, 

i  <  1,  0  <  o;a  <  -TT. 

Therefore,  =  r a  cos  uja,  i"!  =  rQsino;^  is  a  solution  of  and  ri  >  0.  This 
completes  the  proof  of  the  existence  theorem.  □ 

We  just  proved  that  for  sufficiently  large  a  non-parallel  solution  exists.  In  ESI, 
we  showed  that  for  a  <  1  only  parallel  solution  exists.  As  we  reduce  a,  the  branch  of 
non-parallel  solution  cannot  simply  disappear  into  nowhere.  It  must  be  connected 
to  another  branch  of  solution.  The  non-linear  integral  equation  m  is  symmetric 
with  respect  to  ri.  If  (ri,r3)  with  ri  >  0  is  a  solution  of  then  (— ri,r3)  is  also 
a  solution  of  ®.  It  is  natural  to  conjecture  that  as  we  reduce  a  the  non-parallel 
solution  with  ri  >  0  disappears  at  a  =  a*  where  it  is  connected  to  its  counter¬ 
part  (mirror  image)  with  ri  <  0.  This  conjecture  is  a  reasonable  one  and  it  has 
been  confirmed  by  our  extensive  numerical  simulations  (for  example,  see  results  in 
Figure  2).  With  this  conjecture,  a*  is  the  critical  value  for  a  such  that  for  a  >  a* 
non-parallel  solution  exists  while  for  a  <  a*  only  parallel  solution  exists.  In  this 
paper,  we  are  going  to  assume  this  conjecture  is  true  and  we  proceed  to  analyze 
asymptotic  behaviors  of  the  critical  value  a*  as  —  1^0. 

Theorem  5  below  identifies  the  equation  for  a* ,  the  intersection  point  of  the  non¬ 
parallel  solution  (ri  >  0)  with  the  horizontal  line  ri  =  0.  Then  in  the  next  section, 
we  do  asymptotic  analysis. 


Theorem  5.  Suppose  a  non-parallel  solution  (ri  >  OJ  of  Ejl  conneets  to  a  parallel 
solution  (ri  =  0)  at  a*  and  r^.  Then  a*  and  satisfy 

(m3)|r-i=0  =  J's-  (w?)|ri=0  =  ^-  (39) 

a* 

Proof.  We  begin  with  the  fact  that  a,  ri  and  r^  satisfy  system  ®-  Evaluating  the 
third  equation  in  ®  at  ri  =  0,  a*  and  rg,  we  get 

(to3)|^i=o  =  rl- 
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Differentiating  the  probability  density  function  with  respect  to  s  =  a  ri ,  we  have 


-1 


/  miexp[— C/(r 

Js 


i)]di 


exp[— C/(m)]  +  — mi  exp[— C7(m)] 

Zj 


(40) 


=  -(mi)p(m)  +  mip{m). 

At  the  intersection  point,  a(ri)|j.j=o  =  ol*  and  s|ri=o  =  0-  Using  the  definition  to 
calculate  the  derivative,  we  obtain 


dri 


ri—0 


ri  —  0 


lim 

1-1^0  a{ri)ri  —  0 


lim  — — - 

ri^o  a(ri) 


1 

a*' 


(41) 


Here  by  introducing  variable  s  =  a  ri  and  using  the  definition  to  calculate  the  de¬ 
rivative  dri/ ds\r^=Q  directly,  we  avoid  the  assumption  that  the  non-parallel  branch 
of  ri  >  0  and  the  non-parallel  branch  of  ri  <  0  are  connected  smoothly  at  ri  =  0. 

Differentiating  the  first  equation  in  ©  with  respect  to  s  =  ari  and  using 
yields 

^=-(mi)2  +  (m?)  (42) 

Evaluating  it  at  ri  =  0,  using  (I41II  and  (mi)|i.j=o  =  0,  we  arrive  at 

(mi)  |ri=0  =  ~- 

We  remark  that  a*  >  1  is  a  direct  consequence  of  the  above  formula.  □ 


3.  Asymptotic  behavior  of  the  critical  value  a* .  In  order  to  study  the  behav¬ 
ior  of  a* ,  the  critical  value  of  the  intermolecular  dipole-dipole  interaction  strength, 
we  switch  to  the  spherical  coordinate  system  (</>,  0)  where  the  pole  is  on  the  2-axis. 
In  the  spherical  coordinate  system  {(j),6),  we  have 

(mi,  m2,  m3)  =  (sin ())cos0,  sin ())sin0,  cos ())), 

U{(j),9)  =  UMut{(l>,0)  +  UExti4’,&), 

where  the  mutual  interaction  and  the  external  part  of  the  potential  are  given  re¬ 
spectively  by 


Umui  ('(',  0)  =  —ari  sin  (j)  cos  0  —  ar^  cos  (j), 


UExt{4>,&)  =  -pEcoscj)-  ‘t’- 


The  steady  state  probability  density  is 


P(0:^)  =  ^exp[-C7(()),( 


Z  = 


exp[— !/(((),  0)]d0  sin  (j)d(l) 


'  0  ^  —n 
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Lemma  1.  In  the  spherical  coordinate  system  ((/>,  ^);  equation  becomes 
/•i 


liexp 


/-I 


/  77’  I  *  *  N  I  ^0  T— ,2  2 

(/iL  +  a  r^)u  + 


du 


exp 


'-1 


{fj,E  +  a*r^)u  + 


du 


(43) 


—  u^)  exp 

{fiE  +  a*rl)u  +  ‘^E^u^ 

du 

exp 

(/iU  +  a*rl)u  +  ^E^u^ 

du 

1 

a* 


Proof.  In  the  spherical  coordinate  system  {(j),0),  we  introduce  substitution  u  = 
cos  4>.  Then,  we  have 

n-K  /*27r 

i  Vn  ^  ^ 

d9  sin  (j)d(j) 


Z  = 


p~  pZli  j. 

/  /  exp  {piE  +  a^rg)  cos  (p  +  -^E^  cos^  0 

7o  ^0  2  . 


/o  ^0 
27r  J  exp 

pTZ  /*27r 


+  a*r^)u  +  ^E'^u'^ 


du, 


Z{m3)  = 


f! 

/*! 


cos  0  exp  {^lE  P  a* r’l)  cos  (j) +^E‘^  cos^  (j)  d9  sin  (fdtj) 


=  27r  /  wexp 


J-i 

ptz  /*27r 


{fj,E  +  a*r^)u  +  ^E^u^ 


2 

dw, 


Z{ml)  = 


sin^  (j)  cos^  9  exp  {^E  +  a*rl )  cos  (j)+  ^E^  cos^  (j)  d9  sin  (fdcj) 


10  Jo 
/•i 


TT  /  (1  —  exp  (/ri?  +  a*r3)M  + 

J-i  L  2 


2 

drt. 


Substituting  these  into  (P|l.  we  obtain  m- 


□ 


We  consider  the  case  where  oq  and  E  are  fixed  while  {\ao\E  —  fj,)  0. 

Lemma  2.  (A)  When  oq  >  0,  we  have 

lim  -a*(/r)r*(/r)=oo.  (44) 

(B)  When  ap  <  0,  we  have 

lim  _a*(Ai)»'3(M)  =  oo.  (45) 

tJ.^{-ao}E 

Proof.  We  prove  (A)  by  contradiction.  (B)  can  be  proved  in  a  similar  way.  Suppose 
(Ell  is  not  true.  Using  Theorem  3,  we  have  —a*{iP)r^[^)  >  fiE.  The  supposition 
that  (Ell  is  not  true  implies  that  there  is  a  bounded  sequence  of  —a* {ii)r^{ii)  as 
aoE.  Note  that  —  1  <  <  1  is  also  bounded.  As  a  result,  there  is  a 

subsequence  {/in}  satisfying 

lim  /in  =  aoE, 

n — >00 


lim  -a*{fj,n)rl{fin)  =  a, 


(46) 


lim  -r^inn)  =  b, 
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where  both  a  and  b  are  finite:  —1  <  b  <  1  and  a  >  uq  =  .  Taking  the  limit  on 

both  sides  of  yields 


uexp 


/  N  I  2 

(oo  -  a)u  +  —u 


du 


=  -b. 


exp 


/  \  I  2 

(oo  -  a)u  +  —u 


du 


(47) 


/-I 


(1  —  u‘^)  exp  (ao  —  a)u  + 


du 


exp 


'-1 


/  \  I  ^0  2 

(ao  -  a)u  +  —u 


Notice  that  m  and  UHl)  can  be  rewritten  as 


du 


2b 


2b 

a 


(u)  =  -b,  (1  -  u^}  = 

a 

where  the  average  is  taken  with  respect  to  the  probability  density 


p{u)  =  i  exp  (ao  -  a)u  + 


We  consider  the  quantity 

((1  -  u'^)[{ao  -  a)  +  aow]) 
1  '■1 


(48) 


(49) 


Z  = 


exp 


/  \  I  OO  2 

(ao  -  a)u  +  —u 


Oq  2 


=  7T  /  (1  -  W  )[(ao  -  a)  +  aow]  exp  {ao  -  a)u  +  — 


/-I 


du 


i  J  (l-'a^)d|exp  (ao  -  a)-a  +  y I 


du. 


=  -y_^aexp 
From  (HH),  we  have 


(oo  —  a)u  +  y'w^  du  =  2{u). 


((1  -  u‘^)u)  =  —  [((1  -  ■a^)[(ao  -  a)  +  aou])  -  (1  -  M^)(ao  -  a)] 
ao 


=  —  [2(u)  -  (1  -  u^)(ao  -  a)]  =  — 
ao  ao 

Combining  this  result  with  we  obtain 


On  the  other  hand, 


—2b - (ao  —  a) 

a 


2b 

a 


((1  —  u^)(l  +  u))  =  J  {1  —  u^){l  +  u)p{u)du  >  0. 


^=0. 

(50) 

a 

>0. 

(51) 

=  00. 

□ 
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Now  we  introduce  another  Lemma  based  on  Laplace’s  method  before  deriving 
the  asymptotic  behavior  for  a*. 


Lemma  3.  Suppose  f{u)  is  a  quadratic  function  satisfying  /(—I)  =  0  and  g{u)  is 
a  quadratic  function.  As  X  —>■  oo,  we  have 


/(A)  = 


f{u)  exp[— Alt  +  g{u)]du 
I  exp[— Au  +  5(u)](iu 


r\  +  (2/v + +  iHg" + + 3rg']i  + 
/  ^ 


1  +  g'^  +  {g"  +  + ' 

where  all  the  derivatives  are  evaluated  at  —1. 


(52) 


Proof.  Using  the  substitution  m  =  —  1  +  — ,  we  have 

A 


f  V  V 

/  /(-l  +  T)exp[g(-1  + -)]exp(-i;)du 

rm  =  ^ - ■ 

J  exp[g(-l  +  -)]  exp{-v)dv 


Expanding  functions  /  and  we  obtain 


r2X 


/(A)  = 


/"(-I) 


A2 


exp 


A2 


exp{—v)dv 


r2X 


rn 


exp 


'  /"  V 

1  +  - - 

2/' A. 


A2 


exp(—v)dv 


\  ,  fg"  +  g'^\ 


exp{—v)dv 


,  fg''  +  g'^\ 

’+»A+  )  V 


exp(—v)dv 


r2X 


f- 

■'  A 


g"  +  g'^  ,  /"  A 


2/' 


)  A^ 


exp{—v)dv 


^  ,  /^  ,  fg''  +  g'^\ 


10 


exp{—v)dv 


where  all  the  derivatives  are  evaluated  at  —1.  Note  that 

j.2\ 

m  / 

X^ooJo  JQ 

and  the  covergence  is  exponential,  the  result  1521)  follows  immediately.  □ 


p2X  poo 

lim  /  z;*  exp(— z;)dz;  =  /  exp{—v)dv  =  k\ 

X^OoJq  Jq 


We  now  present  the  main  asymptotic  results  for  a*  as  (|ao|U  —  >  0. 
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Theorem  6.  Let  e  =  ^ — r-^.  ^45  e 

\aoE\ 


2  +  -S  + 


-OLqE 


0,  a*  behaves  like 

1 


+ 


(53) 


Proof.  We  present  the  proof  for  the  case  of  oq  >  0.  The  case  of  oq  <  0  can  be 
proved  in  a  similar  way  using  a  lemma  similar  to  Lemma  3. 

For  the  case  of  oq  >  0,  we  derive  the  hrst  term  of  expansion 


*  2 
a  =  - 
e 


(54) 


It  turns  out  that  once  we  obtain  the  hrst  term,  it  helps  us  simplify  signihcantly  the 
derivation  of  subsequent  terms.  The  detailed  derivation  of  the  three  term  asymp¬ 
totic  expansion  is  given  in  the  Appendix. 

Let  A  =  —a*r^  and  g{u)  =  gEu  +  ^E^u^.  Lemma  2  shows  that  as  e  — *■  0  we  have 
A  ^  oo.  Applying  Lemma  3  to  the  hrst  equation  of  (03)  with  f(u)  =  1  -fit,  we  have 


1  = 


T  +2g'-^  +3(g"  -b^'^)—  -b 


A2 


A3 


1  +  +  ig"  +  5^^)^  + ' 


(55) 


Applying  Lemma  3  to  the  second  equation  of  03  with  f{u)  =  1  —  we  arrive  at 


Y  +  i‘^g'  “  1)^  +  [^(5”  +  g'^)  ~  ^g']  yi  + 


A3 


1  +  g'^  +  {g"  +  9'^)^  + ' 


(56) 


Solving  (rg  -b  1)  from  (1^  yields 


1  = 


(1  ~  +  +  l^g'  ~  2(ff"  +  g'^)]  Yi  +  ■ ' 


A2 


1  +  g'j^  +  ig”  +  + ' 


(57) 


‘^g'~  +  3(3"  +  3^^)y3  +  •  •  •  —  (1  —  g')^  +  [33^  —  2(5"  +  5''^)] 


A2 


Equating  (1^  and  03  gives 

j  +  ^g  ^2  '  ^  'A3 

which  in  turn  leads  to 

g'^  =  [g'  —  2(3”  +  3^^)]  +  •  •  ■  ■  (58) 

Taking  derivatives  of  g{u)  =  fiEu  +  -^E’^vf  ^  we  have 

g'(-l)  =  -(aoA^  -  gE)  =  -a^E^e,  g"(-l)  =  aoE^  =  0(1). 

The  right  hand  side  of  (1^  is 

g'-2ig"  +  g'^)  =  -2aoE^  +  Oie). 

The  left  hand  side  of  03  must  match  the  right  hand  side.  Thus  it  follows  that 


A=  -  +  ■ 
e 


(59) 


Equation  03  shows  that  as  A  ^  00  we  have  (—1)-  Therefore,  we  conclude 

that 


*  2 
a  =  - 

£ 


(60) 


922 


HONG  ZHOU,  HONGYUN  WANG  AND  QI  WANG 


□ 

Remark:  This  asymptotic  behavior  of  a*  has  been  confirmed  in  our  numerical 
simulations  (see  next  section). 


Figure  1.  Diagram  of  function  u}{r)  and  function  r(u}).  Along 
Lo{r),  dF/doj  =  0  and  along  r{u}),  dF/dr  =  0.  At  the  intersection 
point  (wq,,  Cq),  the  gradient  of  F  is  zero.  Hence  (wq,,  r^)  is  a  solution 
of  equation 


Figure  2.  The  components  of  the  polarity  vector  ri  and  as 
functions  of  a  with  respect  to  ao  =  1,-1.  Parameter  values  used 
here  are:  fj,  =  0.6,  E  =  1,  a  G  [5,20].  Left  panel:  ao  =  1-  Right 
panel:  ao  =  —1. 


4.  Numerical  results.  We  solve  eq.  0  numerically  to  validate  the  asymptotic 
results  obtained  above.  Figure  2  shows  nonparallel  solutions  for  the  case  of  ao  =  1 
(left  panel)  and  ao  =  —  1  (right  panel).  ri  and  are  shown  as  functions  of  a. 
For  ao  =  1  (left  panel),  the  crtitical  value  is  a*  ~  5.7227  and  the  value  of  is 
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negative.  This  is  consistent  with  the  fact  that  for  ao  >  0  the  third  stationary  point 
of  the  external  potential  is  in  the  same  hemisphere  as  the  negative  direction  of  the 
external  field.  The  parameter  values  used  in  simulations  are  listed  in  the  figure 
caption.  For  ao  =  —1  (right  panel),  the  critical  value  is  a*  «  5.3835  and  the  value 
of  ra  is  positive.  This  is  consistent  with  the  fact  that  for  ao  <  0  the  third  stationary 
point  of  the  external  potential  is  in  the  same  hemisphere  as  the  positive  direction 
of  the  external  field.  Figure  3  shows  ri  and  as  functions  of  a.  For  the  case  of 
ao  =  1  (left  panel),  nonparallel  solution  exists  only  in  a  bounded  interval  of  -E. 
The  lower  limit  of  E  is  explained  by  Theorem  2,  which  says  the  necessary  condition 
for  the  existence  of  non-parallel  solution  is  E  >  This  is  only  a  necessary 

condition,  not  sufficient.  The  actual  lower  critical  value  is  Ki  0.7619.  The 

parameter  values  used  in  simulations  are  given  in  the  figure  caption.  The  upper 
limit  of  E  is  explained  by  Theorem  3,  which  says  for  ao  >  0  another  necessary 
condition  for  the  existence  of  nonparallel  solution  is  ^E  <  a(— ra)  <  a.  Again, 
this  is  only  a  necessary  condition,  not  sufficient.  The  actual  upper  critical  value  is 
^ upper  ~  3.3921.  For  the  case  ap  =  —  1  (right  panel),  nonparallel  solution  exists 
when  E  is  above  a  critical  value,  explained  by  Theorem  2.  For  ao  <  0,  Theorem  3 
does  not  impose  an  additional  condition  for  the  existence  of  nonparallel  solution. 


Figure  3.  The  components  of  the  polarity  vector  ri  and  as 
functions  of  E  with  respect  to  ao  =  1,-1.  Parameter  values  used 
here  are:  /x  =  0.6,  a  =  10,  E  G  [0, 10].  Left  panel:  ao  =  1.  Right 
panel:  ap  =  —1. 

In  order  to  confirm  the  asymptotic  result  (1^  for  a* ,  we  compare  the  asymptotic 
result  from  (1311)  with  the  direct  numerical  solution  of  equation  (1131).  Figure  4 
demonstrates  that  the  asymptotic  expansion  is  accurate  for  both  the  case  of  ap  = 
1  (left  panel)  and  the  case  of  ap  =  —1  (right  panel).  Figure  5  focuses  on  the 
performance  of  the  asymptotic  formula  for  e  near  1  where  the  asymptotic  formula 
is  not  expected  to  yield  accurate  result.  For  the  case  of  ap  =  1  (left  panel)  the 
asymptotic  formula  is  very  good  even  for  e  close  to  1.  The  maximum  relative  error 
is  less  than  3%.  For  the  case  of  ap  =  —1  (right  panel)  the  result  of  asymptotic 
formula  is  still  acceptable  although  for  e  near  1  the  relative  error  is  about  25%. 
Finally,  Figure  6  shows  the  absolute  error  in  asymptotic  formula  if33|i  as  a  function 
of  e.  Here  the  error  is  defined  as  the  absolute  value  of  the  difference  between  the 
numerical  solution  of  equation  63  and  the  asymptotic  result  H53I) .  As  shown  in 
Figure  6,  for  both  the  case  of  ap  =  1  and  the  case  of  ap  =  —1,  the  error  decays  like 
or  faster  as  e  ^  0. 
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Figure  4.  Log-log  plots  of  a*(e)  as  a  function  of  e.  Solid  line 
represents  the  numerical  solution  of  equation  63).  Symbols  repre¬ 
sent  the  asymptotic  result  (P|l.  Left  panel:  ao  =  1;  Right  panel: 
ao  =  —1.  They  demonstrate  that  the  asymptotic  expansion  is  accu¬ 
rate  for  both  the  case  of  ao  >  0  and  the  case  of  ao  <  0.  Parameter 
values  used  here  are:  if  =  1,  /r  G  [0.01,0.99]. 


5.  Concluding  remarks.  In  this  work  we  have  studied  the  equilibria  of  Smolu- 
chowski  equation  for  dilute  solutions  of  dipolar  nematic  polymers  and  magnetic  dis¬ 
persions  under  an  imposed  external  field.  We  establish  the  existence  of  a  critical 
value  of  the  dipole-dipole  intermolecular  strength  above  which  a  class  of  nonparallel 
equilibria  exists  where  the  polarity  vector  is  not  parallel  to  the  direction  of  the  exter¬ 
nal  held.  An  asymptotic  formula  for  the  critical  value  is  established  in  terms  of  the 
other  material  parameters  and  compared  with  direct  numerical  computations.  The 
rigorous  as  well  as  asymptotic  results  for  the  equilibria  of  the  Smoluchowski  equation 
will  provide  valuable  insights  into  the  governing  system  of  pdes  which  is  prototypical 
in  kinetic  theories  for  solution  of  polymers  and  dispersions  HSIIIHIESIEIIII 


Acknowledgements.  The  authors  thank  the  anonymous  referees  for  their  con¬ 
structive  comments  and  suggestions.  H.  Zhou  was  partially  supported  by  the  Naval 
Postgraduate  School  and  AFOSR.  H.  Wang  was  partially  supported  by  the  NSF 
Grant  DMS-0317937.  Q.  Wang  was  partially  supported  by  AFOSR  under  grant 
number  F49550-05- 1-0025,  and  by  the  NSF  Grant  DMS-0204243. 

Appendix.  Asymptotic  expansion  of  a*  with  more  terms 

In  this  appendix,  we  derive  H53II .  Recall  that  A  =  — a*r3  and  g{u)  =  fiEu  +  . 

From  the  previous  analysis  in  Theorem  5,  we  know 

g\-l)  =  -{aoE^  -gLE)  =  0  =  ^oE^  =  0(1). 

This  will  be  very  helpful  in  simplifying  the  expansion  below. 
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Figure  5.  Linear  plots  of  a*(£:)  as  a  function  of  e  for  e  near  1. 
Note  that  the  asymptotic  formula  is  derived  in  the  limit  of  £  ^  0 
and  is  not  supposed  to  be  valid  near  e  =  1.  Solid  line  represents 
the  numerical  solution  of  equation  m-  Symbols  represent  the 
asymptotic  result  15dll .  Left  panel:  oq  =  1;  Right  panel:  ao  =  ~1- 
In  the  left  panel,  it  is  clear  that  the  asymptotic  formula  is  very 
good  for  the  case  of  ao  >  0  even  for  e  close  to  1.  Parameter  values 
used  here  are:  E  =  1,  ^  G  [0.01,  0.99]. 


Figure  6.  Error  in  a{e)  as  a  function  of  e.  Error  is  defined  as  the 
absolute  value  of  the  difference  between  the  numerical  solution  of 
equation  m  and  the  asymptotic  result  itS^.  Left  panel:  ao  =  1; 
Right  panel:  ao  =  —1.  From  the  figure,  we  can  see  that  as  e  ^ 
0,  the  error  decays  like  or  faster.  The  blue  line  and  symbols 
represent  the  error.  The  red  dashed  line  shows  y  =  ce^.  The  red 
line  is  drawn  to  guide  your  eyes  and  to  show  the  behavior  of  error 
as  e  ^  0.  Parameter  values  used  here  are:  E  =  1,  ^  S  [0.01,0.99]. 


Lemma  lA:  Suppose  f{u)  is  a  quadratic  function  satisfying  /(—I)  =  0.  As  A  — *■  oo, 
we  have 


I{\)  = 


f(u)  exp[— Aii  +  g{u)]du 
Ff - 


/  exp[— Alt  +  5(u)]d 


u 


n  +  (2/V  +  ni,  +  [3(5"  +  g'^)f'  +  3/"5U 
+  [12gV'f  +  6g"r]^  +  15g^f^  +  ■■■ 

1  +  +  {g"  +  5^^)^  +  3g'5"^  +  ^9"^^  +  •  •  • 


(61) 
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where  all  the  derivatives  are  evaluated  at  (— !)• 

V 

Proof  of  Lemma  lA:  Using  the  substitution  m  =  —  1  +  — ,  we  have 

A 


f  V  V 

/  /(-l  +  T)exp[g(-1  + -)]exp(-u)(ii 

/(A)  =  ^  ^ 


f  ^ 

J  exp[g(-l  +  -)]  exp{-v)dv 


f  and  g,  we  obtain 


/(A)  = 


A 


A2 


exp 


^  ^  A  2  A2 


exjp{—v)dv 


r2X 


exp 


A2 


exp(—v)dv 


/ 

Jo 


1  + 


Ill 

2/' A 


1  !  9 

'+»A  + 


a  I  „'2 


,  9'9”'v^  ,  9”‘^v^ 

A2  ~A3  + 


e  ^dv 


r2\ 


(9 


"  ^  9'^ 


u2  g'g'' 


exp{—v)dv 


fl'l 


’  1  + 

/  +  2r)j  + 

/  //  1  >2  pH  a 

(’?  +LL 

to 

_ 1 

+ 

(  g'g"  ,  9"  /"  A 
2  2  2/'^ 

3  4 

V  g  V 

W  +  ~1?  + 

e-^dv 


r2X 


\  ,  f9''+9'l  ,  g'g"v^  ,  ^4  ^ 


e  "du 


where  all  the  derivatives  are  evaluated  at  —I.  In  the  above  expansion,  we  have 
used  the  fact  g'  =  0  and  omitted  terms  of  the  order  higher  than  l/A^.  For 

g'^ 

example,  we  omitted  the  term  --nr.  Using 

6  A'^ 

n2X  poo 

lim  /  exp{—v)dv  =  exp(— u)d' 

X^OoJq  Jq 

we  obtain  the  result  m-  □ 


—v)dv  =  k\ 


Theorem  lA:  Let  e  =  As  e  ^  0,  a*  behaves  like 

aoE 


(62) 


Proof:  Applying  Lemma  lA  to  the  first  equation  of  m  with  f{u)  =  l  +  M,  it  follows 
that 


rU  I  = 


I  +  9' ^  +  {9”  +  1a  ’ 


(63) 
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Applying  Lemma  lA  to  the  second  equation  of  m  with  f{u)  =  1  —  one  has 


A 

^  +  (25^  “  +  [^9"  +  9'^)  “  5^]  ^  +  [125'g"  -  65"]  ^  +  9"^-^  +  ' 


,1 


1 


1 


1 


1  +  9\  +  {9"  +  9'  )t9  +  ^9' 9" TT  + 


A3 


A4 


A  '  "  ^A2 

Solving  (rg  +  1)  from  (l(^  yields 
^3  +  1  = 


(1  ~  9')^  +  +  [35^  ~  2(5"  +  5'^)]  —  +  [65"  —  9g'g'']  —  125"^^  +  ■ 


_  A3 

1  +  9'^  +  (5”  +  9’^)-^  +  ‘^9’  9"  +  ‘^9”^-^  + 


Equating  PI  with  PI-  we  have 

^  +  ^^9”  +  9'^)-^  +  125V^^  +  1%"^^  +  ■•• 

=  (1  “  ffO  Y  +  +  [35^  “  2(g"  +  g'^)]  —  +  [6g"  -  9g'g"]  Yg  -  125"^y4  + 


A3 


A4 


which  leads  to 


g'X  =  -2g"  +g'-  2g'^  +  [85"  -  9g' g"\  \  -  12g"^^ 

Recall  £:  =  Substituting 

aoE 

g'i-l)  =  -aoE^e,  /(-I)  =  aoE^ 

into  ra.  we  get 

— cxqE‘^£X  =  — 2o;o-£^^  —  o-qE^e  —  2(q;q_E’^)^£^ 


(64) 


(65) 


(66) 


+3aoA^(l  +  SdoE^e)-^  —  12(aoA^)^  (ffA)^  +  ’  ’  ’  ■ 

We  already  know  eA  =  0(1).  We  seek  a  solution  of  the  form  eA  =  co+cie+C2e^+- 
Substituting  this  expansion  form  in  PI,  we  obtain 

—aoE‘^{co  +  ciE  +  C2s‘^)  =  —2aoE‘^  —  aoE'^e  —  2{aoE‘^)‘^e‘^ 


(67) 


+3aoA^(l  +  3aoA^e)- 


-12(00^2)2  +, 


Co  +  cie 

Equating  coefficients  of  corresponding  powers  of  e  gives 

.A  =  2-i.+  (iaoi^2_3\  2^.... 


(68) 
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Equation  PI  tells  us  that 
(-^3) 


l-i 

1  +  +  ■  ■  ■ 

=  i-i 

1  +  0 

A 

[i+yf +  •••] 

A 

[  VavJ 

=  1- 


2-ie 


=  1  -  - 
2 


1  + 


r'+' 


Combining  this  result  with  PI-  we  conclude  that 


a  =- 


sX 


2-  ie+  ( 

- 

3^ 

SJ 

1  H - 

1-1 

(1  +  3^] 

1  +  ••• 

1 

■  1  /I  ^2  3A  2 

€ 

2--e+(^-aoE  - -j  e  +■•• 

+--J 

□ 
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